Two pairing parameters in superconducting grains 
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Unlike bulk superconductivity, where one energy scale, the energy gap, characterizes pairing 
correlations, we show that in small superconducting grains there exist two different such quantities. 
The first characterizes cumulative properties of the grain, such as the condensation energy, and the 
second single-particle properties. To describe these two energy scales, we define two corresponding 
pairing parameters, and show that although both reduce to the bulk gap for large grains, this occurs 
at different size scales. 
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I. INTRODUCTION 

The question of how superconductivity is affected by 
small sample size, raised already by Anderson in 1959 [1], 
has experienced a recent revival of interest, which started 
with the experimental work of Ralph, Black and Tin- 
kham [2,3]. In these experiments it was shown that small 
superconducting grains, with size much smaller than the 
coherence length, but level spacing d smaller than the 
bulk gap A, have a gap of order 2A in their tunneling 
spectrum, and grains in the regime d > A (to be called 
"ultrasmall" regime) do not. Since the pairing parame- 
ter is the basic quantity in bulk superconductivity, sev- 
eral efforts have been made to define pairing parameters 
which are relevant in the regime of small grains [4-9], 
and reduce to A in the bulk limit. The need for such 
new definitions arises for two apparent reasons: (i) for 
ultrasmall grains, in which d > A, the quantity A has no 
direct physical meaning; (ii) in both ultrasmall grains, 
and small grains for which d ^ A, if the grains are iso- 
lated the appropriate ensemble is the canonical one, in 
which the usual definition of A that is used in the grand 
canonical ensemble [see Ag e. in Eq. (3)] is trivially zero. 
In this paper we show that there is a third, fundamental 
reason for the inadequacy of Ag e. to describe pairing cor- 
relations in small superconducting grains. Unlike the sit- 
uation in bulk superconductors, in which all the various 
superconducting properties can be characterized by one 
energy scale, A, in general there exist two distinct energy 
scales that characterize pairing correlations, and the dif- 
ference between them becomes important particularly for 
small grains. The first such energy scale, which we denote 
by As. p., characterizes single-particle properties, such as 
excitation energies and even-odd effects. The second, de- 
noted as Acu, characterizes cumulative properties, such 
as the condensation energy [see Eqs. (10), (17)], to which 
pairing correlations of all the levels up to lu-q contribute. 
The reason why these two scales are in general distinct 
is that levels far from the Fermi energy, namely those 
with energy |e — i^pj between A and the Debye frequency 
ivu (to be called "not condensed" or "far" levels), make 



a more significant contribution to cumulative properties 
than to single-particle ones. (For a discussion on the role 
of the far levels in superconducting grains and persistent 
currents in normal metal rings see Rcfs. [10,11]). The 
contribution of the far levels to physical properties turns 
out to be proportional to the level spacing; in the large 
grain limit where the level spacing becomes exceedingly 
small, this contribution can thus be neglected. In this 
limit, the single-particle and cumulative properties are 
therefore both determined by the correlations of only the 
"condensed levels", i.e. those within A of Ep, so that the 
two scales Acu and Ag.p. become identical. 

Previous attempts [4,6-9] to define, in terms of pair 
correlation functions, a pairing parameter adequate to 
describe small grains in the canonical ensemble resulted 
in parameters characterizing cumulative properties of the 
grain, but not single particle ones. We shall discuss a 
particular definition for such a cumulative parameter, de- 
noted Acu, which is purposefully chosen such that Acu 
reduces to A in the bulk limit. For small grains, we show 
the correspondence of this definition to cumulative prop- 
erties such as the condensation energy. We then define a 
single-particle parameter Ag.p. and show the correspon- 
dence of this definition to single-particle properties of the 
grain. In the bulk limit both Acu and Ag.p. reduce to A as 
expected, but at a different size scale (see Fig. II). For the 
single-particle properties, we find Ag.p. 3> A in the ultra- 
small regime {d > A), and Ag.p. ~ A for larger grains. In 
contrast, for the cumulative properties we find Acu 3> A 
not only in the ultrasmall regime, but also in an "interme- 
diate regime" [10] in which A > c? > A^/wd. Regarding 
the relation between cumulative and single-particle pa- 
rameters, we find Acu S> Ag.p. in both the ultrasmall and 
the intermediate regimes. Note that both the range of the 
intermediate regime, and the value of the pairing parame- 
ters and the various superconducting properties obtained 
below, depend explicitly on wd- This wn-dependence, as 
well as the very existence of the intermediate regime, are 
direct manifestations of the contribution of the pairing 
correlations of the far levels. Finally, we discuss some 
possible strategies for measuring the single-particle and 
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cumulative quantities in small grains. Throughout the 
paper, except when discussing the even-odd effect, we 
assume for simplicity an even number of electrons in the 
grain. 



II. PAIRING CORRELATIONS IN THE 
GRAND-CANONICAL ENSEMBLE 

We consider the reduced BCS Hamiltonian 

^= H ^jS-<TCja-Ad^'c|+ct_Cj_Cj+ . (1) 

where the second sum (the pairing interaction) is re- 
stricted to energies within of E-p, and +(— ) de- 
note spin up (down). The Hamiltonian (1) is the 
usual BCS Hamiltonian used when discussing supercon- 
ducting grains [9] and its validity is discussed in, e.g., 
Refs. [12,13,9]. For a grain with a given, finite num- 
ber of electrons this Hamiltonian has an exact solution, 
obtained by Richardson [14,15], and independently by 
Gaudin [16]. In the macroscopic limit, where the canoni- 
cal and grand-canonical ensembles produce the same re- 
sults, Richardson's solution reduces to the BCS solution 
[17,18]. As a result of the pairing interaction, the ground 
state of a superconductor is different from the nonin- 
teracting Fermi state. It is a coherent superposition of 
various paired many-body noninteracting states, defined 
as Slater determinant of real one electron wavefunctions. 
The coherence means that the amplitudes for all these 
states in the superposition are real, up to an overall, 
global phase factor. This is true both for the BCS wave 
function in the grand canonical formalism 

|BCS) = n(%+^^i''])|Vac) (2) 
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where bj = Cj-Cj+, and for the exact wave function given 
by Richardson's solution [15]. (In the grand canonical 
formalism, the above coherence relates to noninteracting 
many-body states with a given number of pairs. One can 
add a constant phase between states of different num- 
ber of pairs, which is referred to as the superconducting 
phase) . A particular characteristic of the structure of the 
ground state of a superconductor is that the occupation 
probability of levels above the Fermi energy is nonzero, 
and that of levels below the Fermi energy is smaller than 
unity. A measure of the pairing correlations, which ex- 
ploits both the non Fermi like occupation probability and 
the phase coherence mentioned above, is given by the pair 
amplitude, which in the g.c. ensemble is given by 

Ag.c. = Ad^(cj_Cj+). (3) 

j 

For any many-body BCS-like eigenstate, characterized by 
the set {fja} of the occupation probabilities of the BCS 
quasiparticles, one obtains [19] 



Ag.c.({/,>}) = XdY,niV*{^ - fj+ - fj-) . (4) 

3 

Specifically, for the ground state given in Eq. (2), for 
which fj„ = for all j, we obtain 

Ag.c.(g.s.) = Ad^Ujt;^, (5) 

3 

which gives in the bulk limit Ag.c.(g-s.) = 
a;D/sinh(l/A) = A. The mean occupation of level j 
is given by v'j , and uj + v'^ ~ 1 . Note that the nature of 
Ag e. is cumulative, being a sum of the contributions of 
all levels. As a result of the phase coherence mentioned 
above, all the contributions to the sum in Eq. (5) come 
with the same phase. 
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FIG. 1. The two pairing parameters, Acu (top curve) and 
As. p. (bottom curve), normalized to the value of the bulk gap, 
are sketched as function of iV = u!o /d (proportional to the in- 
verse grain size), for A = 0.12. Ag.p. » A for iV < 250, which 
correspond to Ad — A. Acu > A, As. p. both in the ultrasmall 
regime, and in the intermediate regime, up to A?" ~ 7500 which 
corresponds to X^duo = A^. 



III. CUMULATIVE PAIRING PARAMETER Acu 

We now turn to the definition of the cumulative canoni- 
cal pairing parameter. A natural extension of the concept 
of long range order in the bulk, suggests the following 
definition for a canonical pairing parameter: 

jAcuj' = {Xdf j dridr2F{ri,r2) (6) 

where F{ri,r2) is a function characterizing pairing cor- 
relations, which is given by 

F(ri,r2) = (Vi(ri)Vl(ri)V-(r2)V+(r-2)) - (7) 
(Vt(ri)V+(r2))(Vl(ri)V-(r-2)). 
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We will show that this definition for the cumulative 
canonical pairing parameter is adequate, as it has the 
following properties: (i) it is meaningful for a canonical 
ensemble; (ii) in the bulk limit it is equivalent to the 
grand canonical definition (3) for any given BCS eigen- 
state; (iii) it is related to the condensation energy and 
to the mean occupation of the noninteracting levels not 
only in the bulk limit, but also in the opposite limit of 
ultrasmall grains [see Eqs. (10-11) below]. 

Expanding each of the 'ip operators in the basis of the 
noninteracting single-particle eigenstates, we obtain 

lAeup = {Xdf ^((&l6,) - {clcj+}{clc,-}) . (8) 



Since the terms in each of the brackets are number con- 
serving, they are meaningful in the canonical ensemble, 
which is used below for the evaluation of |AcuP for small 
grains. However, in the macroscopic limit, one can use 
the grand canonical ensemble to evaluate |Acu|^ within 
the BCS approximation. By using the Bogoliubov trans- 
formation and the BCS wavefunctions, we find that for 
any many-body BCS-like eigenstate 

|A,u|2 = {\dfJ2''^vtu*vj{l - /,+ - - /,+ - /,_) 

= |Ag.c.P, (9) 

where the last equality is a result of Eq. (4). The canon- 
ical pairing parameter is equal, in the bulk limit, to the 

grand canonical pairing parameter for any many-body 
state, and therefore both definitions are equivalent in 
this limit. Our definition for A^u differs slightly from 
that given in Ref. [8] . The difference lies in the last term 
in Eq. (8), which results in the exact equivalence to Ag e. 
in Eq. (9). 

We now turn to the opposite limit of ultrasmall grains. 
By examining Richardson's exact solution [14,15], it was 
shown that pairing correlations, however small [i.e. even 
for A -C 1 / In {wu / d)] manifest themselves both in the 
form of the ground state wave function and in a finite 
condensation energy -Bcond [10]. Here we make the con- 
nection between these two effects and the canonical pair- 
ing parameter. In particular, we shall show that to lead- 
ing order in A 
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-^cond 



I A, 



2Xd 



(10) 



where [Acuj^ is evaluated for the exact groimd state 
(this result differs from the known bulk result i?cond = 
A^/(2(i) by the occurrence of A in the denominator, 
which is discussed below). Also, Acu is related to the 
sum of the pairing correlations in all the levels, as are 
reflected in their occupation probabilities, since we shall 
show that 



21n2 



(11) 



XdJ2'^ 



(12) 



with = (b^jbj) and = 1 



The quantity Aq 



which has been defined in analogy to Ag.c, reflects the 
non- Fermi- like mean occupation probability in the many- 
body ground state of the noninteracting single-particle 
levels. 

These results are found as follows: To obtain the value 
of I Acu I ^ ill the ground state we use Richardson's equa- 
tions and expressions for the wavefunctions [15], to write 
the ground state to leading order in A, in terms of the am- 
plitudes of the various noninteracting many-body states 
appearing in it: 



<^g.s.(l,...,iV) = l 

(f>g.,.{l,...,N;^j,k)-- 



\d 



2(efe - €j) 



(13) 



The first line refers to the amplitude that all the levels 
below the Fermi energy are occupied by pairs (noninter- 
acting ground state of the system). The second line is the 
amplitude for the noninteracting many-body state which 
is the same as the noninteracting ground state, except 
for one pair excitation from level j to level k. The am- 
plitudes of all the other many-body noninteracting states 
are zero to first order in A. 

Prom Eq. (13), to zeroth order in A, the two terms in 
Eq. (8) cancel each other. The second term in Eq. (8) has 
no contribution to first order in A. The contribution to 
the first term which is first order in A comes from the fact 
that there is finite amplitude for levels above the Fermi 
energy to be occupied, and is given, using Eq. (13), by 
Xd/2{ej — ei) for each j >,i < (by < (>) we refer to 
levels below (above) the Fermi energy). Since the sum 
in Eq. (8) is unrestricted, we get a factor of two, and to 
leading order in A 
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i<d> 



1 



- Ci 



(14) 



The condensation energy was calculated in Ref. [10]. In 
order to compare it with Acu of Eq. (14) we present here 
its value for a general spectrum, which to leading order 
in A is given by 



-Econd = {Xd)^ ^ 



1 



i<,j> 



2ei - 2ei 



(15) 



This result is obtained directly from Richardson's equa- 
tions, and leads to Eq. (10). 

Evaluating Eq. (14) for equally-spaced spectrum, we 
obtain 



]Acu|^ = 21n2- A^dwD . 



(16) 
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where 

densed levels, and i?cond 



The large magnitude (linear in lud) of this result is due 
to the fact that all the amplitudes in (13) have the same 
phase (which is a consequence of the coherence discussed 
in Sec. 11), so that all the terms in Eq. (14) are added 
with the same sign. 

In Ref.JlO] we have shown that E.^nd ^ Ef^nd + ^col^ 
A^/2(i is the contribution of the con- 
A^cjd is the contribution of 
the far (not-condensed) levels, which can be calculated 
perturbatively. Similarly, |AcuP is a sum of the contri- 
butions of all levels below wd, too, and is related to the 
condensation energy in both the BCS and the perturba- 
tive regimes. Therefore, it is natural to hypothesize that 
|AcuP (similarly to Econd) can to a good approximation 
be written as the sum of a bulk contribution from the 
condensed levels and a perturbative contribution from 
the far levels, i.e. | AcuP ~ 2 In 2\^dujD + A^. This would 
imply that A^u ^ A in both the ultrasmall and inter- 
mediate regimes (see Fig 11). Note, that the expression 
for the condensation energy in Eq. (10) is different from 
the bulk expression of the condensation energy in terms 
of A, namely E^^l\ = A^/2d, by a factor of A in the 
denominator. Since we have shown that in the bulk limit 
Acu = ^g.c. ! we conclude that A^u and the condensation 
energy have a different functional dependence on A. The 
above hypothesis results in 



-^cond — 



lA, 



2/(A)d 



(17) 



where /(A) is a monotonic function of A, which equals A 
for A < 1/lniV and unity for A » l/lnA^ (bulk limit). 

We now turn to derive the relation (11) between Acu 
and the mean occupation probabilities. In the exact 
ground state, is given by 



E 

{/i-/iv|je} 



\'P{fi,-,j,-l 



N 



(18) 



where {/i.-./w | j G} is any configuration of iV levels out 
of the 2N noninteracting single-particle levels, that in- 
cludes the level j. To leading order in A, using Eq. (13), 
we obtain for levels j above Ep 



4(e,- 



(19) 



Using the approximation of constant level spacing, we 
find that 



A^ 

4e, 



(20) 



where ej is measured from the Fermi energy. This im- 
portant result shows that the mean occupation w^(e) is 
proportional to e~^, imlike the usual BCS result, where 
for e ^ A the mean occupation is proportional to e~^. 
Since the occupation probability for a single level is, by 



(20), proportional to the level spacing, this term can 
be neglected in the bulk limit. However, in finite-size 
grains, we find that in both the ultrasmall and interme- 
diate regimes (i.e. A < X^/d^ou), the occupation proba- 
bility for energies of the order of Ud is in fact larger than 
that given by the BCS approximation. 

As a result of Eq. (19), to first order in A we find that 

1. 



and Uj 



I f tv equally 



spaced levels, one therefore obtains from Eq. (20) that 



xdj2', 



(21) 



which yields the relation between Aqcc and A^u given 
in Eq. (11). To summarize: the condensation energy, 
pairing parameter, and mean level occupation all receive 
significant contributions from the weak pairing correla- 
tion of the far ( "not-condensed" ) levels up to the Debye 
frequency. As a result, their magnitude is much larger 
than that given by the BCS approximation, not only in 
the ultrasmall regime, but also in the intermediate regime 
where A^/wd < d < A. 



IV. SINGLE-PARTICLE PAIRING PARAMETER 

As.p. 

Unlike the cumulative properties considered above, 
other superconducting properties, such as the excita- 
tion energy and the even-odd effect (quantified by the 
Matveev-Larkin (ML) parameter, see Ref. [5]), are re- 
lated to the blocking of, say, only one or two levels to 
pairing correlations. Therefore, they do not depend on 
the correlations between all possible pairs of levels in the 
grain, but only on the correlations between these selected 
levels and all the other levels. As a result, their values 
in small grains are much smaller than one would get by 
using the bulk analogy with A^u as the pairing parameter. 
We therefore define 



Ai^,^Xd^{{blbj) + {bib])) 



(22) 



where the sum is over the noninteracting single-particle 
levels and Hs a selected such level. Since one frequently 
deals with the lowest energy levels, we define 



."s.p. 



A! 



s.p. 



(23) 



for i being the level closest to the Fermi energy (for our 
considerations the cases that i is below or above Ep are 
equivalent, and we take i to be below Ep). 

In the bulk limit, using the BCS approximation, we 
obtain for the ground state 



Aip=2Xdu,v,Y,'' 
j 

and specifically 



(24) 
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As.p. = Xd^UjVj = A. (25) 

i 

We now turn to the ultrasmall regime, and evaluate As. p. 
in the ground state to second order in A. Using Eq. (13) 
we obtain 




and for the equally spaced spectrum 

As.p. = Ad+A^dln^ . (27) 

Let us now consider the excitation energy, say i?cxc, of 
the first excited state of an (even) superconducting grain. 
This state can be described by having the two levels near- 
est to Ep singly occupied with probability unity, thus 
breaking one pair, and the other N — 1 pairs occupying 
the remaining 2N — 2 levels according to Richardson's ex- 
act solution [15]. The excitation energy has two different 
contributions, (i) The kinetic energy cost d of occupying 
a level of higher energy with probability unity, and (ii) 
a pairing energy cost, which is given by Ag.p.. Note that 
the latter has (iia) a diagonal part, given by Xd, which 
is related to the excess energy of two electrons occupy- 
ing the same spatial noninteracting eigenstate, and (iib) 
an off-diagonal part due to the blocking of the two singly 
occupied levels to pairing correlations. Adding these con- 
tributions gives 

Eexc = d+ A,.p. « d -h Ad + A^d In ^ , (28) 

where the last expression is valid for a grain with equally 
spaced spectrum. Similarly, the ML parameter [5] is 
given by 

Aml = ^A,.p. . (29) 

Note that the pairing contribution to the energy cost 

of the first excitation equals 2As.p. in the bulk limit, 
whereas it equals Ag.p. in an ultrasmall grain, and the 
same factor of two appears in the ML parameter. This 
reflects the fact that in the bulk limit each level is cor- 
related with all the other levels, while in the ultrasmall 
grain, the dominant correlations are between levels on 
different sides of the Fermi energy. 

V. MEASUREMENT STRATEGIES FOR Acu AND 



A possible measurement of the cumulative correlations 
of a superconducting grain was discussed in Ref. [10]. It 
was shown that the condensation energy can be obtained 
from specific heat or spin magnetization measurements, 
where an explicit calculation was done for the latter. 



In distinction with the case of bulk superconductivity, 
where the single-particle properties are easiest to measure 
through the lowest excitations of the system, in small 
grains the energy of the first excited state is given by 
Eq. (28), in which Ag.p. is only a small correction to the 
level spacing. Matveev and Larkin suggested to measure 
their pairing parameter by the parity-induced alternation 
of Coulomb blockade peak spacings in small grains. Here 
we suggest the possibility that for an ensemble of small, 
weakly-coupled grains, Ag.p. could be measured through 
their spin magnetization as a function of magnetic field 
at zero temperature. We assume that the coiipling be- 
tween the grains is weak enough such that the equilib- 
rium properties of the system can be approximated by 
summing over the individual grains, but strong enough 
so tunneling between grains in a large enough ensemble 
(of say 50 grains) occurs within the time of the experi- 
ment. 

For a single-grain, the spin magnetization shows steps 
at values that correspond to the Zecman energy being 
equal the energy required to break a pair, of which the 
first is at i^exc = d + Ag.p. [by Eq. (28)], which is dom- 
inated by the large single-grain level spacing d. In or- 
der to avoid having to worry about the latter, we con- 
sider an ensemble of weakly connected superconducting 
grains, which would have an effective joint level spac- 
ing dons that is much smaller than d if the ensemble is 
large, dens d (we neglect the charging energy. Note, 
that the charging energy was found to be much smaller 
[20,21] than what naive estimates give. Since the tun- 
neling between the grains is weak, the relevant ensemble 
is the canonical one). Then the ground state of the un- 
connected system is given by each of the grains having 
an even number of electrons, being in its ground state. 
The first excited state of the system is given by moving 
one electron between two grains, so that afterwards the 
first has a single-occupied level just below Ep, and the 
second has a singly-occupied level just above Ep. For a 
system of "normal" noninteracting grains the energy of 
such an excitation is of the order of the single-particle 
level spacing of the whole ensemble of grains, dens- How- 
ever, in a system of ultrasmall superconducting grains, 
since two singly occupied levels are created in two dif- 
ferent grains, the energy of such an excitation would be 
given by Ag.p. (which is 3> dens)- Since the first step in 
the spin magnetization is at the energy of the first excited 
state of the system, in such a measurement we predict a 
gap of value Ag.p.. (In a similar measurement in grains 
with d < A, the measured gap would be given by A). 
Thus, Ag.p. can be interpreted as the superconducting 
gap as measured by single particle properties. In princi- 
ple, performing such a measurement as a function of the 
grain size would monitor the change of Ag.p., from being 
much larger than A for ultrasmall grains to equaling A 
for large ones. Note by that reducing the grain size one 
increases this gap Ag.p.. 

Finally, we would like to note that the logarithmic 
dependence of the correlation (second) term of Ag.p. in 
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Eq. (27) is manifested in the interaction correction to the 
ensemble averaged magnetic response of small metallic 
grains, when considered within the BCS model [11]. 



VI. SUMMARY 

We have shown that various superconducting proper- 
tics can be classified to two groups, one containing those 
properties which are single particle in nature, and the 
other containing those properties which are cumulative, 
a result of the summed contributions of many levels. Un- 
like the case in bulk superconductivity, where both these 
properties are characterized by one energy parameter A, 
in general, and in particular in small grains, two different 
energy parameters, Ag.p. for the former group and A^u 
for the latter one, characterize superconducting proper- 
ties. 

In bulk superconductivity the relevant contribution to 
the various superconducting properties comes from the 
"condensed" levels, within A of Ep. However, in small 
grains the contribution of all levels up to wd is significant, 
and this results in the existence of these two different pa- 
rameters, as well as in the fact that Acu ^ Ag.p. ^ A in 
the ultrasmall regime and Acu ^ A also in the interme- 
diate regime. Experimental possibilities to measure the 
two pairing parameters were discussed. 
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